Phase-dependent light propagation in atomic vapors 
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Light propagation in an atomic medium whose coupled energy levels form a ^-configuration 
exhibits a critical dependence on the input conditions. Depending on the relative phase of the input 
light fields, the response of the medium can be dramatically modified and switch from opaque to 
semi-transparent. These different types of behaviour are caused by the formation of coherences 
due to interference in the atomic excitations. Alkali-earth atoms with zero nuclear spin are ideal 
candidates for observing these phenomena which could offer new perspectives in control techniques 
in quantum electronics. 
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I. INTRODUCTION 

Experimental evidence has demonstrated that the non- 
linear optical properties of laser-driven atomic gases ex- 
hibit counter-intuitive features with promising applica- 
tions. A peculiarity of these media is the possibility to 
manipulate their internal and external degrees of free- 
dom with a high degree of control. Recently the control 
of the internal dynamics in an atomic vapor by means of 
electromagnetically induced transparency (EIT) [l[ was 
demonstrated for the generation of four- wave mixing dy- 
namics 0] and of controlled quantum pulses of light 
Zeeman coherence has also been used to induce phase de- 
pendent amplification without inversion in Samarium va- 
pors [5[ and in HeNe mixtures [6|]. In another experiment, 
the interplay of internal and external degrees of freedom 
in an ultracold atomic gas by means of recoil-induced res- 
onances Q was used to achieve waveguiding of light Q . 
From this perspective, it is important to identify further 
possible control parameters on the atomic dynamics for 
the manipulation of the non-linear optical response of the 
medium. 

Recent studies have been focusing on the dynamics 
of light interacting with atoms featuring coupled en ergy 
levels in a so-called 'closed-loop' configuration 0, fiof . 
In this configuration a set of atomic states is (quasi-) 
resonantly coupled by laser fields so that each state is 
connected to any other via two different paths of co- 
herent photon-scattering. As a consequence, the rela- 
tive phase between the transitions critically influences 
dynamics @ and steady states [l(| EH, EH- Applica- 
tions of closed-loop configurations to nonlinear optics 
have featured double-A systems where two stable or 
metastable states are -each- coupled to two common ex- 
cited states. A rich variety of nonlinear optical phenom- 
ena has been predicte d fill. [13, Ihi a nd experimentally 
observed @, H S, El Efl U3, El SSl- In Hi, in par- 
ticular, it has been shown experimentally that the prop- 
erties of closed-loop configurations can be used to cor- 
relate electromagnetic fields with carrier frequency dif- 
ferences beyond the GHz regime. Moreover, coherent 
control based on the relative phase in closed-loop con- 



figuration has been proposed in the context of quantum 
information processing [2tj . 




FIG. 1: Electronic transitions of the ^-configuration. Each 
transition ji) — > \j) is resonantly driven by a laser field at fre- 
quency Uij. Here, \g) is the ground state, |1) and |2) the inter- 
mediate states, which decay into the ground state at rates 7i 9 
and 72 S , respectively, and |e) the excited state, which decays 
with rates 7 e i and 7 e 2 into the corresponding intermediate 
states. Each pair of levels is coupled by two paths of exci- 
tation, hence the dynamics depends critically on the relative 
phase between the paths. The coherent dynamics of the 0- 
configuration is equivalent to that of the double-A scheme, 
whereas the radiative instability of the atomic levels differs. 



In this work we investigate the phase-dependent dy- 
namics of light propagation in a medium of atoms whose 
energy levels are driven in a closed-loop configuration, 
denoted by the <0> (diamond) scheme and depicted in 
Fig. [TJ This configuration consists of four driven tran- 
sitions where one ground state is coupled in a V-type 
structure to two intermediate states, which are in turn 
coupled to a common excited state in a A-type structure. 
It can be encountered, for instance, in (suitably driven) 
isotopes of alkali-earth atoms with zero nuclear spin [2l| . 
Although the coherent dynamics of <0 schemes is equiva- 
lent to that of double-A systems [9(, the steady states of 
the two systems exhibit important differences due to the 
different relaxation processes [lj], E2] • 

The dynamics of light propagation in a medium of <0>- 
atoms is studied by integrating numerically the Maxwell- 
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Bloch equations. We find that, depending on the in- 
put field parameters, the polarization along the medium 
can be drastically modified. The propagation dynam- 
ics may exhibit two metastable values of the relative 
phase, namely the values and 7r, corresponding to a 
semi-transparent and to an opaque medium, respectively. 
For different values of the initial phase, light propagation 
along the medium tends to one of these two values, de- 
pending on the input values of the driving amplitudes. 
These two types of the medium response are supported 
by the formation of atomic coherences leading to a mini- 
mization of dissipation by depicting the population of one 
or more atomic states. This phase dependent behavior, 
selected at the input by the operator, offers promising 
perspectives in control techniques in quantum electron- 
ics. 

The article is organized as follows. In Sec.|TT]the model 
is introduced and discussed. In Sec. lIIII the results for the 
dynamics of light propagation, solved numerically from 
the equations reported in Sec. Ill Dt are reported and dis- 
cussed in some parameter regimes. Conclusions and out- 
looks are reported in Sec. IIVI The appendices present 
in detail equations and calculations at the basis of the 
model derived in Sec. [Til 



II. THE MODEL 

We consider a classical field propagating in a dilute 
atomic gas along the positive z-direction. The field is 
composed of four optical frequencies u\ g , v-i g , v e \ and 
v e 2, its complex amplitude is a function of time t and 
position z of the form 

E(z,t) = ±Y,£ ij (z,t)e ij e- i ^ t - k « z+ ' l >^ z ' t » +c.c, 

where fey denotes the wave vector and the polariza- 
tion of the frequency component i/y. The input field 
enters the medium at z = 0, and the effect of coupling 
to the medium is accounted for in the z dependence of 
the amplitude £ij(z,t) and phase <fiij(z,t) whose varia- 
tions in position and time are slow with respect to the 
wavelengths Ay = 2-K/kij and the oscillation periods 
T = 2ir/vij, respectively. The atomic gas is very di- 
lute and we can assume that the atoms interact with the 
fields individually. In particular, each field component 
at frequency i/y drives (quasi-) resonantly the electronic 
transition \i) — > \j) of the atoms in the medium, such 
that the atomic levels are coupled in a <)-shaped config- 
uration. 

The relevant atomic transitions and the coupling due 
to the lasers are displayed in Fig. [TJ The ground state 
\g) is coupled to the intermediate states |1), \2) at ener- 
gies tuv\ , hbJi by transitions with dipolc moments di ff = 
(l|d|g) and d 2g = (2|d|g), respectively. The intermedi- 
ate states decay back into the ground state at decay rates 
7i s and 72^. The intermediate states are also coupled to 



the excited state |e) at an energy fku e with respect to the 
ground state \g), by the dipole transitions d e i = (e|d|l), 
d e2 = (e|d|2). The excited state |e) decays into states 
|1) and 1 2) at rates 7 e i and r y e 2, respectively. A similar 
configuration of levels can be found in isotopes of alkali 
atoms with zero nuclear spin [2lj . 

The light fields propagating through the dilute atomic 
sample will induce a macroscopic polarization in the 
atoms. This polarization will depend on intensities and 
phases of the light fields. The polarization, in turn, will 
affect absorption and refraction of the light fields, alter- 
ing their propagation. Below we introduce the equations 
for field propagation and the corresponding atomic dy- 
namics. 



A. Equations for field propagation 

We denote by P(z,i) the macroscopic polarization in- 
duced in the atomic gas 

P(z,t) = nTr{da(z,t)} (2) 

where d is the dipole operator, n is the density of the 
medium, which we assume to be zero for z < and uni- 
form for z > 0, and a(z, t) is the atomic density matrix at 
time t and position z, which has been obtained by trac- 
ing out the other external degrees of freedom. Details of 
the underlying assumptions at the basis of Eq. ^ are 
discussed in Appendix A. 

We decompose the polarization P(z, t) into slowly- and 
fast-varying components, namely 

P{z,t) = l^2v ij ( z ,t)e ij e- i ^ t - k » z +^^ +c.c, 

whereby the complex amplitudes Vij and the phases <pij 
vary slowly as a function of position and time. We con- 
sider the parameter regime where the driving fields are 
sufficiently weak so that the generation of higher-order 
harmonics can be neglected. By comparing Eqs. @ 
and ([3]), the amplitudes Vij can be expressed in terms 
of the elements of the atomic density matrix a, 

Vij = 2nV lj a lj e l( ^ t - k ^ z+x ^ ) (4) 

where <7y = (i|<r|j). We have expressed the dipole mo- 
ments in direction of the electric field polarization as 
• dji — T>ije~ lSi:i , thereby separating the complex am- 
plitudes Vij into modulus and phase. Here, the term 
is real, Oij are the dipole phases (#y = —dji), and 

Xij{z ) t) = (f>ij(z ) t)-O i j (5) 

is the sum of the slowly- varying field phases <^>y (z, t) and 
the dipole phases 0y. 

Using definitions fTJ) and and applying a coarse- 
grained description in time and space, the Maxwell equa- 
tions simplify to a set of propagation equations for each 
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of the slowly- varying components of the laser and polar- 
ization fields 1221 



d£jj 
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"" 1 -Re{P ij (e kh cf> kl )}, (7) 



2 £q c £., 



which are defined for z > 0. Here, each amplitude £ij and 
phase (f>ij is coupled via the corresponding polarization 
Vij to all other field amplitudes and phases. 

We rescale the propagation equations using the dimen- 
sionless length and time 



£ = KlgZ, T = CKigt. 

Here K\ g is the absorption coefficient 



(8) 



, 1 vigVlg 



such that determines the characteristic length at 

which light driving the transition \g) — ► |2) penetrates 
a medium with density n. We denote the dimensionless 
field amplitudes by 



where 



Q,ij(z,t) =T>ij£ij(z,t)/h 



(9) 



(10) 



is the real valued Rabi frequency for the transition \i) — > 
\j). In this notation the propagation Eqs. ^ and (J7|) 
reduce to the form 



dQjj , 8Q, 



8t 



where 
Pij it, T ) 



o-y-exp 



CK lg 



K1 S 



(11) 

(12) 
(13) 



denotes the atomic density matrix elements in a rotated 
reference frame. In the remainder of this paper we con- 
sider laser field geometries where |1) and |2) are states 
of the same hyperfine multiplet so that v\ g ~ v^g and 



B. Atomic dynamics 

The time evolution of the density matrix u{z, t) for the 
atomic internal degrees of freedom at position z > is 
governed by the master equation 



& = — [H{z,t),(x] +£<t. 



(14) 



where z is a classical variable. Equation (fT4|) is obtained 
by tracing out the degrees of freedom of momentum and 
of position in the transverse plane, in the limit in which 
the medium is homogeneously broadened and the atoms 
are sufficiently hot and dilute such that their external 
degrees of freedom can be treated classically. Details of 
the assumptions at the basis of Eq. (fT4|) are reported in 
Appendix A. Here the Hamiltonian 

H{ z ,t) = H-Ij'XjI ( 15 ) 

- J Yl (V jg (z,t) e-^-^+x^V \j)(g\ 

3=1,2 

+fl ej (z,t) e -*("^*-*«i*+w(*.*)) \ e )(j\+B..c. 

describes the coherent dynamics of the internal degrees 
of freedom, and it depends on z through the (real- valued) 
Rabi frequency ttij(z, t) given in Eq. dTUJ) , and through 
the field and dipole phases, Eq. 

The states |1), |2) and |e) are unstable and decay radia- 
tively with rates 7i g , j2g and ~f e — j e \ +7e2, respectively. 
The relaxation processes are described by 

La = E ^^kMm-\m°-°\m\) 

3=1,2 

+ W^Mem - \e)(e\a ~ <r\e)(e\) , (16) 

3=1,2 

where the recoil due to spontaneous emission is neglected 
since the motion is treated classically. In the remainder of 
this paper we assume a symmetrical decay of the excited 
level, 7 e i = 7 e2 = 7 e /2. 

We note that the transitions \g) — > \j) (j — 1,2) are 



saturated when 
transitions \j) — > \ 



For later convenience, we introduce 



Tjg. Correspondingly, the upper 

ej > le + log- 



> 

are saturated when f2 P1 > 



Gej = 



l+7e/7; 



(17) 



jg 



which explicitly shows the scalings of the upper field am- 
plitudes with the corresponding decay rates. 



C. The relative phase 

In so-called closed-loop configurations, like the <^> 
scheme, transitions between each pair of electronic levels 
are characterized by -at least- two excitation paths, in- 
volving different intermediate atomic levels [lfl, [3 ■ In 
the <) scheme the relative phase between these excita- 
tion paths critically determines the solution of the master 
equation, and hence the atomic response during propaga- 
tion. The role of the relative phase in the atomic response 
is better unveiled by moving to a suitable reference frame 
for the atomic evolution, which is defined when all am- 
plitudes £ij are nonzero. 
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We denote by p the density matrix in this reference 
frame, obeying the master equation 



1 

P = T 
in 



H,p 



Cp. 



In this reference frame the Hamiltonian |T 
formed to H, El 



(18) 
is trans- 



H 



HA e \e)(e 
\ (fiei|e)(l 
+fii fl |l)(sl 



ftAi|l)(l| + fiA 2 |2>(2| 
V Sl e2 e* 

«2 S |2)( 9 | 



(19) 



j9(*,t) | e ><2j 
H.c), 



with the detunings 



Al = U)\ - V\g, 

A e = UJ e - U e i - fig. 



(20) 
(21) 
(22) 



The Hamiltonian (|19p exhibits an explicit dependence on 
the phase 



Q(z, t) — Av t - Afc z + Ax(z, t). 



(23) 



where 



Av = v el + v lg - v 2g - v e 2, (24) 
Afc = fc el + fc iff - fc 2ff - fc e2 , (25) 

&X(Z, t) = Xel (z, t) + Xlg{*, *) - Xe2(z, t) ~ X2g(z, t), 

(26) 

with Xij as defined in Eq. (0. 

The four-photon detuning Av results in a time- 
dependent phase, the wave- vector mismatch Afc in a posi- 
tion dependent phase, and Ay f z, t ) comprises the relative 
dipole and field phases. In [IJ, [I^] it has been discussed 
how Q(z, t) affects the dynamics and steady state of the 
atom. The latter exists for Av = and in the remainder 
of this article we assume 

Av = 0, Afc = 0, 

i.e., the atoms are driven at four-photon resonance and 
by copropagating laser fields, such that the wave vector 
mismatch is negligible. Hence, the phase 

@{z,t) = A X (z,t) 

depends solely on the relative dipole phase, which is 
constant, and on the relative phase of the propagating 
fields, which evolves according to the coupled Eqs. (TlT]) 
and ([121) . 



D. Propagation of the field amplitudes and phases 

Having introduced the basic assumptions and defini- 
tions, we now report the equations for the propagation of 
the field amplitudes and phases in the <^>-medium, which 



are numerically solved in Sec. Mil We relate the elements 
of the density matrix p in the new reference frame with 
the elements p io from eq. H31) by p gX = p g i, p g2 = p g2 , 
Pel = Pel, and 

Pe2 = Pe2exp (~i0) . 

The propagation equations for the light fields in the new 
reference frame can then be obtained from Eqs. (| 1 1 [) - ((T2|) 
and take the form 



")9 



Re{p jg } 



for j = 1,2 and 



dQei j s x 
-j—- = -\m{p e i} , 

8<f> e i _ Re{p e i} 



/el 



eg, 



c2 



f = -lm{ P e 2 e M } 



e2 _ Re{ P e2e lB } 



(27) 
(28) 

(29) 
(30) 
(31) 
(32) 



d? Qe2 

where we have introduced the variable 
C = i + r. 

These equations describe the evolution of field ampli- 
tudes and phases as a function of the atomic density ma- 
trix elements pij. In turn, the values of pij depend on 
the field amplitudes and the relative phase O according 
to Eqs. (fTS|) and ([!§)) . The propagation dynamics now 
can be investigated by solving the coupled Eqs. (TTS)) and 
(f2"T)) - ([3"2"|) . The optical Bloch equations for the density 
matrix p are presented in Appendix B. 

In general, the density matrix elements entering 
Eqs. (|2"7| - (f3"2")) are time dependent, i.e., p = p(r). In this 
paper we consider the case of sufficiently long laser pulses, 
such that the characteristic time of change of amplitude 
and phase of the fields and the interaction time between 
light and atoms exceed the time scale in which the atom 
reaches the internal steady state. In this regime, we 
can neglect transient effects, and the density matrix ele- 
ments entering Eq. (f2"?| - ({3"2"|) are the stationary solutions 
of Eq. (TIB)) satisfying dp/dt = 0. This assumption allows 
us to neglect the time derivative in Eqs. (|2"7l) - ([32"|) , hence 
taking £' w £. 

The numerical study of the solutions of Eqs. J27p - (13"2")) 

presented in this paper is restricted to certain parame- 
ter regimes that single out the role played by the phase 
and the radiative decay processes in the dynamics. In 
particular, we consider the situation where each atomic 
transition is driven at resonance, namely 

A, = 0, 



5 



for i — l,2,e. Moreover, we restrict ourselves to the 
regime where the fields are initially driving the corre- 
sponding transitions at saturation. This latter assump- 
tion is important to guarantee a finite occupation of the 
excited state |e), and thus to highlight the dependence of 
the dynamics on the relative phase 9. 

During propagation, it may occur that one of the field 
amplitudes vanishes in just one point of the propagation 
variable When this happens, the relative phase 9 is 
not defined and its value has to be reset manually by im- 
posing continuity of the trajectory, when integrating the 
field equations in amplitude and phase (see Eqs. ([2"T]) - 
(f32| ). The correctness of this procedure has been checked 
by comparing the results with those obtained by integrat- 
ing the field equations for the real and imaginary parts 
of the complex field amplitudes. 



III. LIGHT PROPAGATION IN THE 0-MEDIUM 

In this section we summarize some peculiar proper- 
ties of the <0-level scheme, which have been extensively 
discussed in [T2|. These properties provide an impor- 
tant insight into the propagation dynamics, which we 
study by solving numerically the Maxwell-Bloch Equa- 
tions, Eqs. (fT8|) and (f27 |) - ([32|) . in the regime where the in- 
put fields couple resonantly and saturate the correspond- 
ing electronic transitions, as described in Sec. Ill Dl 



A. Symmetries of the level scheme 

Before entering the detailed discussion of the numerical 
results, it is instructive to review some basic properties 
of the <C>-level scheme, which significantly affect its re- 
sponse to light-propagation. Special symmetries of this 
configuration are encountered when the laser amplitudes, 
resonantly driving the upper (lower) transitions, are ini- 
tially equal, namely when 



Gel = Ge2 = Ge 
Gig = G2g = Gg- 



(33) 
(34) 



In this regime, the Hamiltonian Q19[) substantially sim- 
plifies for some values of the phase. In particular, for 

9 = 0,tt 

(modulus 2tt) the dynamics can be mapped to those of 
well-known three-level schemes [12J. Insight is gained 
by studying Hamiltonian (|19p in a convenient orthogonal 
basis set of atomic states. 

For 9 — ir, Hamiltonian ([IT)]) describes coherent cou- 
pling within the orthogonal subspaces {\g), l^ii)} and 
{|e),|*r 2 >}» = [|l)±|2)]A/2. These sub- 

spaces are decoupled: state l*^) ^ s decoupled from |e) 
and \^i 2 ) from \g) by destructive interference. Sponta- 
neous decay eventually pumps the atom into the sub- 
space {\g), l 1 !^)}, see Fig.^a). Hence, in the stationary 



regime the excited state is depleted and the atomic levels 
which scatter light can be mapped onto a V-level scheme. 

For = and under condition (|33|) . ([34]) . the Hamil- 
tonian (|19[) can be written in two equivalent ways: It 
can describe coherent scattering among the orthogonal 
levels {|1), l^^g), |2)}, forming a A-configuration, while 
state I't'eg) i s decoupled, or alternatively, it can de- 
scribe coherent scattering among the orthogonal levels 
{\g}, 1*12)5 l e )}i forming a S-level scheme, while state 
^12) is decoupled. Here, |*^) are symmetric and anti- 
symmetric superpositions of the states |e) and \g). How- 
ever, if spontaneous decay is included, the two schemes 
are not equivalent. The relaxation processes select one 
configuration over the other depending on the stability of 
state \^eg)y which decays at a rate j e , with respect to the 
stability of state l*^), which decays at a rate ji g + 72^. 
It is then important to introduce the parameter 



7e 



7lff + l2g 



(35) 



which is the ratio between the decay rates of the two 
decoupled states, or, equivalently, the ratio between the 
decay of the excited and intermediate states. Hence, for 
a <C 1 (i.e. the excited state is longer lived than the 
intermediate ones), l*^) * s essentially empty, and the 
effective dynamics can be mapped to a S-level scheme, 
see Fig. EJb). For a ^> 1 instead (i.e. the intermedi- 
ate state is longer lived than the excited one), \^eg) ^ s 
empty, and the effective dynamics can be mapped to a 
A- level scheme, see Fig. (He). Hence, if the ratio a is 
sufficiently different from unity, the dynamics of the <0>- 
scheme can be mapped to three-level schemes and leads 
to coherent population trapping (CPT) [l[ in the sta- 
tionary state. In the A case (a 3> 1), a large coher- 
ence between the intermediate states is observed as re- 
ported in the transient dynamics of pulse propagation in 
a medium of <0>-atoms [15j . Here, for some parameter 
regimes one can observe population inversion at steady 
state on the transition \g) — + |1), |2) 12]. In the S case 
(a <C 1), a macroscopic coherence between ground and 
excited states is created. For some parameter regimes 
one can observe population inversion at steady state on 
the transition |1), |2) ->• |e) [Hj]. 

These properties have important consequences for the 
propagation dynamics. We note that for 9 = 0, tt 
the components of the polarizations Re(pi g ), Re(p2g), 
Re(/?i e ), Re(p2e exp(iO)) vanish. This means that the 
field phases remain constant upon propagation in agree- 
ment with Eqs. J35J, §5U§ and Hence, if at the input 



e(f = o) = o, 7T, 



then 



99 







(36) 



(37) 



and the relative phase remains constant during propaga- 
tion along the medium. We recall that for 9 = it we 
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FIG. 2: Level schemes onto which the <0 configuration is mapped, if the lower field amplitudes, as well as the upper field 
amplitudes, are equal. For O = it, the <0 system imitates a V-configuration (a). Mapping to the E-configuration, (b), is 
obtained for = and when |e) is metastable. Similarly, the mapping to the A-configuration (c) is obtained for = and 
metastable intermediate states. 



observe a V-type dynamics (from now on denoted as de- 
structive interference) and for 9 = metastable CPT on 
a S or A-scheme (from now on denoted as constructive in- 
terference). Hence, from these simple considerations we 
expect that for different values of the input phase and re- 
laxation rates, energy will be dissipated at very different 
rates along the medium. 



B. Destructive interference in the atomic 
excitations 

For 0(0) = 7r the atoms are perfectly decoupled from 
the upper fields independently of their intensity and the 
upper state is empty as described in IIII Al Destructive 
interference makes the polarizations of the transitions be- 
tween the intermediate and the upper states as well as 
the population of the excited state to vanish identically, 
i.e., p e i = Pe2 = Pee = [II] . Correspondingly, the dy- 
namics of light propagation of the lower fields is expected 
to be that encountered in a medium of V-atoms. 

Figure [3Ja) displays the propagation dynamics along 
the medium for O(0) = ir and equal initial field ampli- 
tudes, Gij{0) = Go- Here, one sees that the upper fields 
propagate through the medium as if it were transparent, 
keeping a constant value. The amplitudes of the lower 
fields display identical decays. Figure 0(b) presents the 
corresponding populations of the energy levels along the 
medium. The energy level |e) remains depleted while 
the intermediate states |1) and |2) maintain the same 
population as a function of £ corresponding to the fact 
that the lower fields decay identically along the medium. 
The value of ground and intermediate state populations 
is the saturation value of the corresponding dipole tran- 
sition until about £ ~ 200 when the lower fields Gj g (£) do 
not saturate the transition any longer. After this pene- 
tration length only the ground state is appreciably occu- 
pied. Note that these dynamics are independent of the 
upper field amplitudes, as they remain decoupled from 
the atoms. 



We can find an analytic expression for the dynamics 
shown in Fig. [3J and for the propagation length of the 
lower fields by solving the propagation equations (|19[) 
and ([H])- (32]) for 6(0) = = n. Setting Gej = G e and 
Qj g = Q gi we obtain the equations for the dimensionless 
amplitudes 



8Ge 



l + 4<? 2 



= 0. 



(38) 
(39) 



Here the right hand side in Eq. (|39|) vanishes since 
lm{pij(Q = 7r)} = 0. Therefore, the relative phase © 
and the upper field amplitudes Ge are constant along the 
medium and the medium is transparent for the upper 
fields. Equation ([35)1 is the equation for an electric field 
propagating in a medium of resonant dipoles so that the 
lower field amplitudes G g decay during propagation at 
a rate that depends only on the value of G g itself. In 
the case of large input intensities (see Fig. [3]) a simple 
equation for G g (0 is obtained [23| 



GA0 = ^G 2 g {Q)-m 



(40) 



allowing for an estimate of the penetration depth 
(2Gq(0)) of the lower fields in the medium. 



Constructive interference in the atomic 
excitations 



As discussed in section IIII Al for = the response 
of the system is similar to that of a 3 or of a A level 
scheme, depending on the ratio of the decay rates a in 
Eq. ([35| . Atomic coherences between either the interme- 
diate states or the ground and excited state may form, 
and correspondingly the imaginary part of the polariza- 
tions may become very small, thus reducing dissipation. 

Figure 2] displays the propagation dynamics along the 
medium for different values of the ratio a for 0(0) = 0. 
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FIG. 3: (color online) Propagation of the phase (a), field am- 
plitudes (b) and the corresponding atomic state populations 
(c) for an input phase of 0(0) = n. Here, Q e is the rescaled 
amplitude as in Eq. (JT7J. The phase is constant and the up- 
per fields propagate unperturbed through the medium, while 
the excited state |e) remains depleted. The behavior is inde- 
pendent of Q e and 7 e and thus of a. 
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FIG. 4: (color online) Field amplitudes as a function of the 
propagation length £ for input parameters = 0, Q e — Q g = 
1 and different ratios between the decay rates: long lived 
intermediate states, a = 100, in (a), balanced decay rates, 
a = 1, in (b), and long lived excited state, a — 0.01, shown 
in (c). Correspondingly, the phase O = remains constant 
along the medium (not shown). The rate of dissipation is 
critically determined by a and is slower for a sufficiently larger 
or smaller than unity. 



For a ^> 1 and a -C 1 the amplitudes decay slowly as 
a function of £, as expected from the formation of EIT- 
coherences. Figure [S] shows light propagation for the case 
when the initial conditions of the fields give rise to popu- 
lation inversion at steady state due to metastable CPT. 
We find that population inversion is maintained until 
£ ~ 100 along the absorbing medium, but it gradually 
decreases, since the atomic coherences that are support- 
ing CPT are not stable. 

In the simulations of Figures 2] and the lower (upper) 
field amplitudes remain equal during propagation. If we 
assume that Q e j — Q e and Qj g = Q g for all relevant £, then 



the propagation equations for the amplitudes reduce to 

(41) 



dGe 



Qe Gn a 

9 (1 + 2a), 



Do 
Do 



[G 2 g + a + 2a 2 + Ql (1 + a)] , (42) 



with 



D = Gt (1 + a) + (1 + 4G 2 g ) a {G 2 + a + 2a 2 ) 

+G't [a (2 + 3a) + & (1 + 3a + 2a 2 )] . (43) 



Equations (|4Tj) and (|42j) describe the dissipative propa- 
gation of the field amplitudes, and exhibit a nonlinear 
dependence on the amplitudes and the ratio a of the de- 
cay constants. Here, one can see that for different values 
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e = o, g e = i, g g 



10, and q = 10. Population inversion 



between the intermediate levels and the ground level is found 
along the medium until £ ~ 100. 



of a the absorption lengths can vary by orders of mag- 
nitude. Limiting cases are found for a — > 0, i.e. when 
the excited state is stable, and for a — ► oo, i.e., when the 
intermediate states are stable. In these cases, the right 
hand sides of Eqs. (flTj) and (|4*2")1 vanish, damping is ab- 
sent, and light propagates through the medium as if it 
were transparent 24 1 . 



D. 



Stability of interference under amplitude 
fluctuations 



In the cases discussed so far, the input phase is a con- 
stant of the propagation, and the amplitudes of the upper 
fields, as well as the amplitudes of the lower fields remain 
equal along the medium. We now address the question 
of stability of these configurations against phase and am- 
plitude fluctuations. 



Numerical investigations show that at 8 = ir the V- 
level dynamics is robust against phase and amplitude 
fluctuations, from which we infer that this is a stable 
configuration. It should be remarked that the overall 
behavior is transient in that the medium dissipates the 
lower fields until (well inside the medium) the atoms are 
all found in the ground state. 

The case 6 = is more peculiar. In the cases dis- 
cussed in Sec. IIII CI energy is exchanged between upper 
and lower fields until the lower field amplitudes go below 
saturation. Then, the upper fields decouple as the pop- 
ulation of the intermediate states becomes negligible. In 
order to study the long term dynamics of propagation, we 
now focus on the regime where the lower transitions are 
driven well above saturation and where we may expect 
different length scales for the propagation of the upper 
and lower fields. 

Figure [5] displays propagation when the lower transi- 
tions are driven well above saturation for different values 
of a. The dynamics observed in the a = 1 case separates 
the regimes corresponding to a S-like response for a«l, 
and a A- like response for a 3> 1. For the separating case 
of a = 1 in Fig. [SJc) and (d) we find that the damp- 
ing of the lower fields below saturation is accompanied 
by a drop of the population from intermediate to ground 
states. For a = 100, shown in Fig. [^a) and (b) propaga- 
tion is characterized by EIT-like coherences between the 
intermediate states, which are established through the 
medium by the action of the lower fields. These coher- 
ences increase the penetration depth of the lower fields 
in the medium and allow the upper fields to propagate 
quasi undamped. This is consistent with the behaviour 
discussed in Sec. IIII CI However, for a — 0.01 in Fig.[Bfd) 
and (e) we observe a clear deviation from the symmetric 
decay of the upper field amplitudes at long propagation 
length. 

We now focus on this case, which exhibits novel fea- 
tures with respect to the cases studied so far. In Fig.[^e) 
one observes that the upper field amplitudes, Q e \ and Q e 2 1 
initially equal to each other, undergo a transient behavior 
where they become different: energy is transferred from 
one field to the other, such that the amplitude of one in- 
creases while the other gradually vanishes. This behavior 
is accompanied by a depletion of the excited state, while 
the intermediate states continue to be equally populated. 
At the same time of the vanishing of one of the upper field 
amplitudes, the phase jumps from to 7r and energy 
is redistributed between the two upper fields till they 
reach almost the same value. After this transient, the 
field amplitudes of the excited states remain at a con- 
stant value across the medium. Correspondingly, during 
and after this transient, the excited state population in 
Fig. decreases until it reaches zero. This remark- 
able behavior hints to an instability of the phase value 
= 0, which seems to be triggered here by numerical 
fluctuations of the values of the upper field amplitudes. 
Such conjecture is supported by the numerical analysis 
shown in Fig. where we have introduced fluctuations 
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FIG. 6: (color online) Propagation of fields (upper row) and corresponding atomic states populations (lower row) as a function 
of the propagation length and for different values of the excited state decay rate. The initial conditions are 0(0) = 0, Gg(0) = 10 
and t? e (0) = 1. Here, (a) and (b) correspond to the case a = 100, (c) and (d) to the case a = 1, and (e) and (f) to the case 
a — 0.01. In this latter case, the upper field amplitudes, C? e i and Q e 2, become different: energy is transferred from one field 
to the other, such that the amplitude of one increases while the other gradually vanishes. At this point, the phase jumps 
to the value 7r, and energy is redistributed between the upper fields till they reach the same value. This behavior hints to an 
instability of the phase value O = 0, which seems to be triggered by numerical fluctuations of the values of the upper field 
amplitudes. 



between the initial values of the upper field amplitudes. 
As the initial discrepancy increases, the splitting of the 
upper field amplitudes appears at earlier locations in the 
medium although the behavior of the lower fields is unaf- 
fected. More detailed investigations on populations and 
phases show that the imbalance between the upper field 
amplitudes induces a depletion of the excited state un- 
til the vanishing of one of the upper amplitudes forces a 
phase jump to the value = -k and the upper fields de- 
couple from the atom. After the phase jump, the upper 
field amplitudes tend to recover an equal value, but they 
decouple from the atoms once the lower field amplitudes 
have vanished. 

An explanation of the phase jump from O = to = tt 
is the tendency of the system to minimize the rate of 
dissipation in a way reminiscent of what is observed in 
Four- Wave Mixing experiments where interference effects 
are gene rated in order to minimize spontaneous emis- 
sion [26j . We also note that with increasing values of a, 
the splitting of the upper fields for the same initial dif- 
ference in their amplitudes is progressively delayed inside 
the medium and eventually disappears. 



E. Generic phase at the input fields 

Having identified and investigated two special values 
of the input phase, we now address the question, how 
the phase evolves starting from a generic input value, 
and correspondingly, how light propagates and is dissi- 
pated along the medium. We restrict to the configuration 



with initially equal lower field amplitudes and equal up- 
per field amplitudes (|33[) and (|34|) , and vary the input 
phase 8(0) — tt/2 in steps of 7r/4. 

Figures [5] and [5] display the amplitude and the rela- 
tive phase of the propagating fields for different values 
of the excited state decay rate: a = 100 and a = 0.01. 
Although the lower field amplitudes are clearly damped 
for all values of a, the mechanism of radiation dissipa- 
tion depends on a and on the initial strengths of the field 
amplitudes. This can be associated with a particular evo- 
lution of the phase along the medium, which in the cases 
displayed in Fig. [5] reaches the stable value = 0, and 
in the cases displayed in Fig. [5] tends first to the value 
7r before eventually reaching = 0. The choice between 
these behaviors depends on the input amplitudes of the 
fields. We now discuss these two behaviors in detail. 

In Fig. [5] all atomic transitions are driven at satura- 
tion, and the saturation of the upper transitions is larger 
or equal to that of the lower transitions. We observe 
that the relative phase of the fields tends to the zero 
value. Before this value is reached, radiation is damped 
at a fast rate. Once = 0, the rate of damping of the 
lower field amplitudes changes abruptly to a lower value. 
This sudden change occurs at a propagation length deter- 
mined by the typical absorption length of the fast decay- 
ing transition. The system simulates a V-configuration, 
thereby switching to an EIT-like response. A similar kind 
of behavior is also observed in a medium of the Double- A 
atoms where EIT-coherences are established between the 
two stable states [ll|. In the § configuration the coher- 
ences and interferences are transient 161. For a slower 
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decay of the excited state, however, the system can also 
switch to a S-dynamics and exhibit a transient coherence 
between ground and excited states. A manifestation of 
this phenomenon is population inversion between the ex- 
cited and the intermediate states along the medium in 
Fig. Of). 

In Fig.[H]the lower transitions are driven well above sat- 
uration, and the corresponding saturation parameter is 
larger than the saturation parameter of the upper fields. 
Here, during a transient regime the phase slowly tends 
to 6 = 7r. Nonetheless, the tendency of the medium for 
long propagation lengths is to eventually decouple upper 
fields and atoms, i.e. to switch to a V-dynamics. The 
onset of this dynamics depends critically on the value of 
Q g , which must well saturate the lower transitions with 
respect to Q e in order to populate the intermediate states 
on a time scale shorter than their decay rate, but long 



enough for incoherent decay of the upper state to take 
place. This behavior is in agreement with Sec. IIII D[ 
showing that when the lower transitions are driven well 
above saturation the value = ir is the only stable phase 
under amplitude and phase fluctuations. 



F. Four- wave mixing 

So far, we have considered input fields with equal lower 
and upper field amplitudes. We now discuss propaga- 
tion when one field is initially very weak while the other 
three transitions are driven at or above saturation, and 
study how the dynamics of energy redistribution among 
the fields depends on the input parameters and on the 
stability of the excited state. 

Figures [TO] and QT] display the field propagation when 
the upper field amplitude Q e 2 is very small and the phase 
is initially set to the value 0(0) = 7r/2. In both figures 
we have assumed the excited state to decay slower than 
the intermediate states, but one may also observe ampli- 
fication of the weak field under different conditions. In 
Fig. [TO] the three input fields drive the respective tran- 
sitions well above saturation. Here, amplification of the 
weak field is accompanied by the asymptotic approach- 
ing of the phase to = it. The field Ge2 is amplified 
until the upper state is depleted because of interference 
between the upper fields. From this point further the 
phase = 7r is stable and the lower fields dissipate, until 
they drop below saturation. The jump of the phase to 
the value is an artifact due to all atoms being in the 
ground state. 

In Fig. [TT] the three input fields Q\ g ,Q2 g and Q e \ are 
just saturating the respective transitions. Here, ampli- 
fication of the field Q e 2 is accompanied by a transient 
stabilization of the phase at = 7r. This is accompanied 
by a fast decrease of the lower field amplitude Q g2 , until, 
when Q g 2 vanishes, the phase falls to = 0. After this 
point the behavior changes and G g 2 first increases slightly 
and then decays slowly as a function of £ in a way similar 
to G g i, while the upper field amplitudes remain constant. 
The final configuration supports a coherence between the 
excited and the ground state, and indeed for = and 
this value of a the dynamics can be mapped to a S-level 
scheme. In particular, due to destructive interference, 
the fields are only weakly coupled to the transitions and 
the medium is semitransparent. This is also supported 
by Fig.fTlTc) where one sees that the population is redis- 
tributed between the ground and the excited state while 
the intermediate states are depleted. In this regime the 
medium is characterized by population inversion between 
the excited and the intermediate states. 



IV. DISCUSSION AND CONCLUSIONS 

We have investigated numerically light propagation in 
a medium of atoms whose electronic levels are resonantly 
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FIG. 8: (color online) Propagation of relative phase (a), field amplitudes (b), and populations of the atomic levels (c), for 
various input phases O(0) = 0, ±tt/4, ±7t/2, ±37r/4, it. The upper fields are saturated at least as deeply as the lower fields. 
Curves (a-c) are obtained for input parameters a = 100, and Q g = Q e — 10. Curves (d-f) have input parameters a = 0.01, 
Q g = 5 and Q e = 10. During propagation the phase tends to O = 0. Once this phase is reached, the rate of energy dissipation 
along the medium changes abruptly to a significantly lower level. In the case(d-f) this change is accompanied by the establishing 
of population inversion on the transition |1), 2) — > |e), see (f), due to the formation of EIT-coherences. 



driven by lasers in a configuration. Propagation is crit- 
ically affected by the initial parameters of the input fields 
and shows the tendency to reach configurations which 
minimize dissipation. An important role is played by 
the relative phase 9 between the fields. It exhibits two 
fixed points, = and = it, whose stability during 
propagation depends on the field amplitudes and on the 
ratio a between the rates of dissipation of excited and 
intermediate states. A generic input phase evolves, in 
general, to one of these values, again depending on the 
input amplitudes and a. 

These two metastable phase values are associated with 
two different types of atomic coherences. The response of 
the medium, corresponding to the phase = 0, is char- 
acterized by the formation of atomic coherences typical 
of EIT-media. Similar behaviors have been observed for 
instance in [l5|, [l6| and are analogous to the response 



predicted for light propagation in double- A media [ljj . 

The response of the medium for the phase = n is 
supported by a different type of interference which leads 
to a depletion of the upper state and to a complete decou- 
pling of the upper fields from the atom. For this value of 
the phase, the medium acts like a V- level configuration. 
We note that this value of the phase appears to be the 
preferred value for the <0 medium if the lower transitions 
are driven well above saturation. This behavior is novel 
to our knowledge and it is reminiscent of the phenomenon 
of suppression of spontaneous emission observed in four- 
wave mixing studies in atomic gases [26j . 

In general, the system exhibits a rich dynamics and 
several novel features due to atomic coherence which offer 
new perspectives in control techniques in quantum elec- 
tronics. These could be studied in atomic gases where the 
ground state has no hyperfine multiplet, like e.g. alkali- 
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FIG. 9: (color online) Propagation of relative phase (a), field amplitudes (b), and populations of the atomic levels (c), for various 
input phases O(0) = 0, ±vr/4, ±7r/2, ±37t/4, n. Curves (a-c) are obtained for input parameters a = 5, and Q g — 10, Q e = 1, 
and curves (d-f) for a = 0.5 and the same initial amplitudes. Here, the lower transitions are driven well above saturation and 
during a transient regime the phase tends to the value 7r, while the upper fields decouple. The transition to the value O = at 
large values of £ is an artifact, since for these lengths the lower fields are very weak and the atoms are essentially in the ground 
state. 



earth isotopes which are currently investigated for atomic 
clocks 21] . 

In the future we will extend our analysis to the case 
in which the transitions are not resonantly driven and 
we will address the asymptotic behavior of the system 
following the lines of recent works [27], [28| . 
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APPENDIX A: MACROSCOPIC POLARIZATION 
IN THE SEMICLASSICAL LIMIT FOR THE 
ATOMIC MOTION 

We consider the dynamics of the density matrix g 
of the atomic internal and external degrees of freedom, 
where the center-of-mass degrees of freedom are treated 
as classical variables. Hence, the position x and momen- 
tum p are parameters, distributed according the func- 
tion u>(x, p) which we assume to be stationary, with 
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FIG. 10: (color online) (a) Propagation of relative phase (a), 
field amplitudes (b), and atomic populations (c) for input 
parameters O = 7r/2, a = 0.1, Q g \ = Q g 2 = Qie = 10, and 
Qie = 0.1. Energy is exchanged between fields t?2e and Q g \ 
and also between Qi e and Q g 2 until the excited state decouples 
and the upper fields propagate freely. The jump of the phase 
to the value is an artifact due to all atoms being in the 
ground state. 
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FIG. 11: (color online) (a) Propagation of relative phase (a), 
field amplitudes (b), and atomic populations (c) for input 
parameters O = 7r/2, a — 0.01, G g i = Q g 2 = Gie = 1, and 
t/2e = 0.01. The amplitude t/2e while the phase approaches 
= 7r. At the same time Qi e gradually vanishes at which 
point the phase changes to = and the fields propagate 
almost undamped. 



J dxdpw(x, p) = N and N is the number of atoms. The 
spatial density of atoms n(x) is found from w(x., p) ac- 
cording to Jdpw(x, p) = n(x). In this work we assume 
uniform density, namely 

n(x) = n 

with n constant. The master equation for the density 
matrix g, at the point (x, p) in phase space has the form 

Q=-[H{x,p;t),Q]+£Q (Al) 

where the Hamiltonian H(n,p;t) contains the coherent 
dynamics of the atoms driven by the classical field, 

H{x,p;t) = £-=+H(z,t) (A2) 



and H(z,t) is defined in Eq. (fT5]) . The Liouvillian C 
in Eq. (fll))) describes the relaxation processes, which we 
consider here to be purely radiative. The corresponding 
macroscopic polarization has the form 



P(x, t) = / dp^(x, p) Tr{dp(x, p)} (A3) 



Assuming that the atomic gas has been Doppler cooled, 
so that line broadening is homogeneous, the kinetic en- 
ergy can be neglected in evaluating the atomic response 
to the light. By integrating over p and x, y we hence ob- 
tain Eq. (fbi]). whereby a(z) — J dpdxdyw(x, p)er(x, p), 
and polarization as in Eq. @. 
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APPENDIX B: OPTICAL BLOCH EQUATIONS where p i3 = p*^, p gg = 1 - p ee - p n - p 2 2, and we have 
IN THE PHASE- REFERENCE FRAME taken j el = 7 e2 = 7 e /2. 



We consider Master Eq. (fl8|) in the reference frame of 
the phase. With the notation /5 e2 — p e 2 exp(— «8), the 
corresponding optical Bloch equations are given by 

Pee = i~^(Ple ~ Pel) + i^(/5 2e (Bl) 

~Pe2) - lePee 

&le , \ , .^lg i 

Pn = (Pie -Pel) — (Pgi - Pig) 

+yPee-7i ff Pn (B2) 

.^2e/~ ~ \ , .^2g , 

P22 = -l—(P2e-Pe2)+l — (Pg2-p2g) 

+ y Pee ~ l2gP22 (B3) 
Pel = (l (A X -A e )- ^^P 1 ) Pel (B4) 



. ^le / \ , .^2e s& .^ln 

(Pll -Pee) +« — e P21-«— / 



p e2 = ( i (A 2 - A,) _ 7e + 729 ) p e2 ( 155 ) 



, . ^2e, \ , .^le _iO .^2q 

+«— (P22 "Pee) +1 — e Pl 2 -« — e p e g 

Pi 9 = - («Ai + + i-J-Peg ( B6 ) 

.^2g . .^lg , n 

— P12 + I — (Pgg - Pll) 

p 2g = - (*A 2 + 23a) p 2g + i^e^Vg (B7) 

.^lg . • ^2g, \ 

— P21 +« — (Pgg -P22) 



,i, 2 = ( i (A 2 -A 1 )-^^)p 12 +^e^ e 2 

-i%Pg 2 -^ l9 -^ Ple (B8) 



7e\ . .^le 



— 6 P2g - ^Pel - l^fe U Pe2- (B9) 
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